Preferential interaction coefficient for nucleic acids and other cylindrical poly-ions 
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The thermodynamics of nucleic acid processes is heavily affected by the electric double-layer 
of micro-ions around the polyions. We focus here on the Coulombic contribution to the salt- 
polyelectrolyte preferential interaction (Donnan) coefficient and we report extremely accurate ana- 
lytical expressions valid in the range of low salt concentration (when polyion radius is smaller than 
the Debye length). The analysis is performed at Poisson-Boltzmann level, in cylindrical geometry, 
with emphasis on highly charged poly-ions (beyond "counter- ion condensation"). The results hold 
for any electrolyte of the form z_:z+. We also obtain a remarkably accurate expression for the 
electric potential in the vicinity of the poly-ion. 



Coulombic interactions between salt and poly-anions 
play a key role in the equilibrium and kinetics of nucleic 
acid processes [![. A convenient quantity quantifying 
such interactions and allowing for the analysis and inter- 
pretation of their thermodynamics consequences, is the 
so called preferential interaction coefficient. Several def- 
initions have been proposed and their interrelation stud- 
ied, see e.g. In the present work, they are defined 
as the integrated deficit (with respect to bulk conditions) 
of co-ions concentration around a rod-like poly-ion. Our 
goal is to provide analytical expressions describing the 
effect of salt concentration and poly-ion structural pa- 
rameters on the preferential interaction coefficient, for a 
broad class of asymmetric electrolytes. For symmetric 
electrolytes, it will be shown that our formulas improve 
upon existing analytical results. For other asymmetries, 
they seem to have no counterpart in the literature. Our 
analysis holds for highly (i.e. beyond counter-ion conden- 
sation and uniformly charged cylindrical poly-ions, 
and is explicitly limited to the low salt regime (i.e. when 
the poly-ion radius a is smaller than the Debye length 
1/k). These conditions are most relevant for RNA or 
DNA in their single, double, or triple strand forms. 

As in several previous approaches 0, H, H EU, we 
adopt the mean-field framework of Poisson-Boltzmann 
equation, in a homogeneous dielectric background of per- 
mittivity e. The same starting point has proven relevant 
for related structural physical chemistry studies of nu- 
cleic acids ll|. In a electrolyte, the dimensionless 
electrostatic potential = eip/kT (with e > the ele- 
mentary charge and kT thermal energy) then obeys the 
equation [l2j 
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where r is the radial distance to the rod axis. The va- 
lencies z + and z_ of salt ions are both taken positive. 
Denoting derivative with a prime, the boundary condi- 
tions read r<j)'(r) — 2£ > at the polyion radius (r = a) 
and ^ — > for r — > oo. The latter condition expresses 
the infinite dilution of poly-ion limit and ensures that 
the whole system is electrically neutral, since it (indi- 
rectly) implies that rip' — > for r — > oo. We consider 
a negatively charged poly-anion for which <p < and 
the line charge density reads A = —e^/ts < 0, where 
£b = e 2 /(ekT) denotes the Bjerrum length (0.71 nm in 
water at room temperature). Finally, the Debye length is 
defined from the bulk ionic densities and through 
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interaction coefficient is defined as [7|, [8|. 
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while its cationic counterpart follows from electro- 
neutrality. This quantity -which provides a measure of 
the Donnan effect [3]- can be expressed in closed form 
as a function of the electrostatic potential, see Appendix 
lAl As can be seen in (|A3I) and (|A4[> . T depends expo- 
nentially on the surface potential <po, so that deriving a 
precise analytical expression is a challenging task. Fur- 
thermore, we are interested here in the limit na < 1 
(including the regime na <^ 1) which is analytically more 
difficult than the opposite high salt situation where to 
leading order, the charged rod behaves as an infinite 
plane, and curvature corrections can be perturbatively 
included dEGllIl- 

We will proceed in two steps. Focusing first on the 
surface potential <po = 4>{a), we make use of recent re- 
sults [18| that have been obtained from a m app in g of Eq. 
(T]) onto a Painleve type III problem [H, [2CJ, HU • The 
exact expressions thereby derived only hold for 1:1, 1:2 
and 2:1 electrolytes, but may be written in a way that is 
electrolyte independent. This remarkable feature is spe- 
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z+/z 111/10 1/3 1/2 



1 



10 



c 



-2.51 -1.94 -1.763 -1.502 -1.301 -1.21 -1.06 



TABLE I: Values of C appearing in Eq. Q as a function 
of electrolyte asymmetries. For z+/z~ = 1, 1/2 and 2, C is 
known analytically from the results of [l8| ] . The corresponding 
values are recalled in Appendix[B] For other values of z+ / z- , 
C has been determined numerically, see in particular Fig. HJ] 
of Appendix [Bl 



cific to the short distance behaviour of <j> and has been 
overlooked so far, since not only short distance but also 
large distance properties have been studied [HI]. We are 
then led to conjecture that the corresponding expression 
holds for any binary electrolyte z_:z + , and we explicitly 
check the relevance of our assumption on several specific 
examples. 

Technical details are deferred to the appendices. It is 
in particular concluded in Appendix [B] that the surface 
potential may be written 
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Expression Q is valid for na < 1 and z + £ > 1 [in fact 
z + £ > 1 + 0(1/| log na |)]. These conditions are easily 
fulfilled for nucleic acids. The "constant" C appearing 
in ([3]) depends smoothly on the ratio z + /z- but is oth- 
erwise salt and charge independent. We report in Table 
U its values for several electrolyte asymmetries. The de- 
crease (in absolute value) of C when z + /z_ increases is 
a signature of more efficient (non-linear) screening with 
counter-ions of higher valencies. 



From Eq. ([3]) and the results of Appendix [B] our ap- 
proximation for r takes a simple form 
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This expression is tested in Figures Q] and [5] against the 
"true" numerical results that serve as a benchmark. In 
Fig. Q] which corresponds to a monovalent salt (or more 
generally a z:z electrolyte), we also show the prediction 
of Ref. [9J], which is, to our knowledge, the most accurate 
existing formula for a 1:1 salt. For the technical reasons 
discussed in Appendix [Bl and that are evidenced in Fig- 
ure El our expression improves that of Shkel, Tsodikov 
and Record [9(, particularly at lower salt content. For 
1:2 and 2:1 salts, we expect Eq. ([5]) to be also accurate, 
since it is based on exact expansions. The situation of 
other salt asymmetries is more conjectural (see Appendix 
IB]), but Eq. (O is nevertheless in remarkable agreement 



with the full solution of Eq. (Q}, see Fig. [2] To be spe- 
cific, in both Figures [1] and [H the relative accuracy of 
our approximation is better than 0.2% for na = 10~ 2 
(for both ss and ds RNA parameters). At na = 0.1, the 
accuracy is on the order of 1%. 

As illustrated in Fig. [3l approximation (j4|) assumes 
that z+£ > 1. The corresponding expression for Y there- 
fore breaks down when £ is too low. More general expres- 
sions, still for Ka < 1, may be found in appendix [Cl The 
inset of Fig. [3J offers an illustration and shows that the 
limitations of approximation ([4]) may be circumvented at 
little cost, providing a quasi-exact value for T. Moreover, 
it is shown in this appendix that for z + £ = 1, Jl reads 
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On the other hand, Eq. J3J still holds. The corresponding 
r is shown in Fig. |U 

We provide in Appendix [C] a general expression of the 
short scale (i.e valid up to nr ~ 1) radial dependence 
of the electric potential, see Eq. (|C1[) . The bare charge 
should not be too low [more precisely, one must have 
£ > £c with £ c given by Eq. (|C5p]. and pt -which encodes 
the dependence on £- follows from solving Eq. (|C2[) . In 
general, the corresponding solution should be found nu- 
merically. However, one can show a) that Jl vanishes 
for £ = £ c , b) that Ji takes the value © when z + £ = 1 
and c) that Ji is given by Q when z + £ exceeds unity 
by a small and salt dependent amount. In practice, for 
DNA and RNA, we have £ > 2 and Eq. g]) provides 
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FIG. 1: Preferential interaction coefficient for a 1:1 salt. The 
main graph corresponds to ss-RNA with reduced line charge 
£ = 2.2 while the inset is for ds-RNA (£ = 5). The circles 
correspond to the value of ^ following from the numerical 
solution of Eq. ([1]). The prediction of Eq. ([5]) with Jl given 
by (j4]) and C ~ —1.502, shown with the continuous curve, is 
compared to that of Ref. 0, shown with the dashed line. As 
in all other figures, the opposite of F is displayed, to consider 
a positive quantity. 
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FIG. 2: Same as FigureQ]for a 1:3 and a 3:1 electrolyte. From 
Table HI we have C ~ —1.21 in the 1:3 case and conversely 
C ~ —1.94 in the 3:1 case. The symbols correspond to the 
numerical solution of Eq. ([lj and the continuous curves show 
the results of Eq. (JSJ) with again Jl given by @ • 
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FIG. 3: Preferential interaction coefficient for a 1:1 salt (hence 
C ~ —1.502) and sa = 10 -2 . The circles show the numerical 
solution of PB theory fT} , the continuous curve is for ([5]) with 
Q and the dashed line is the prediction of Ref . Q . Although 
approximation Q breaks down at low £, the inset shows that 
Jl following from the solution of Eq. (|C2|) gives through © a 
P (continuous curve), that is in excellent agreement with the 
"exact one", shown with circles as in the main graph. 



excellent results whenever na < 0.1. To illustrate this, 
we compare in Figure [5] the potential following from the 
analytical expression (|C1[) to its numerical counterpart. 
We do not display 1:1, 1:2 and 2:1 results since in these 
cases, Eq. (|C1[) is obtained from an exact expansion and 
fully captures the r-dependence of the potential. For the 
asymmetry 1:3, Fig. [5] shows that the relatively simple 
form (|C1[) is very reliable. A similar agreement has been 
found for all couples z_:z + sampled, with the trend that 
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FIG. 4: Same as Fig. [T]for £ = 1 and z+/z- = 1. The same 
quantities are shown: our prediction for F [Eqs. ([5]) and (|6]) 
with C ~ —1.502] is compared to that of Ref. The inset 
shows — z+T/z- for a 1:2 salt such as MgCi2 where C takes the 
value -1.301. Circles : numerical data; curve : our prediction. 



the validity of (|C1|) extends to larger distances as z + / z_ 
is decreased. In this respect, the agreement shown in Fig. 
[5]for which is quite high (3), is one of the "worst" 
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FIG. 5: Opposite of the electric potential versus radial dis- 
tance in a 1:3 electrolyte with na — 10~ 2 . The continuous 
curve shows the prediction of Eq. (|C1[) with Jl given by (|4)| ; 
the circles show the numerical solution of Eq. (TTJ. The po- 
tential for £ = 2.2 is shown in the main graph on a log-linear 
scale, and on a linear scale in the lower inset. The upper inset 
is for £ = 5. 

Conclusion. The poly-ion ion preferential interaction 
coefficient V describes the exclusion of co-ions in the 
vicinity of a polyelectrolyte in an aqueous solution. We 
have obtained an accurate expression for Y in the regime 
of low salt (na < 1). The present results are particu- 
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larly relevant for highly charged poly- ions ( z+£ > 1, that 
is beyond the classical Manning threshold [22|]), but are 
somewhat more general and hold in the range £ c < £ < 1 , 
where £ stands for the line charge per Bjerrum length and 
£ c is a salt dependent threshold, given by Eq. (|C5[) . Our 
formulae have been shown to hold for arbitrary mixed 
salts of the form z_:z + (magnesium chloride, cobalt hex- 
amine etc). They have been derived from exact expan- 
sions valid in 1:1,1:2 and 2:1 cases, from which a more 
general conjecture has been inferred. The validity of this 
conjecture, backed up by analytical arguments, has been 
extensively tested for various values of z + /z_, poly-ion 
charge and salt content. These tests have provided the 
numerical value of the constant C reported in Table [H 
which only depends of the ratio z + /z_. As a byprod- 
uct of our analysis, we have obtained a very accurate 
expression for the electric potential in the vicinity of the 
charged rod (r < 

It should be emphasized that the validity of our 
mean-field description relying on the non-linear Poisson- 
Boltzmann equation depends on the valency of counter- 
ions (z + ), and to a lesser extent to the value of z_ fl2l.l23l|. 
For the 1:1 case in a solvent like water at room temper- 
ature, micro-ionic correlations can be neglected up to a 
salt concentration of 0.1M [8|. For z + > 2 or in sol- 
vents of lower dielectric permittivity, they play a more 
important role. Our results however provide mean-field 
benchmarks from analytical expressions, from which the 
effects of correlations may be assessed in cases where they 
cannot be ignored (see e.g. § for a detailed discussion). 
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Combining both relations with adequate weights, in order 
to suppress the integral over counter-ion (+) density, we 
have 
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where (f>o = <fi(a) is the surface potential. Equation (|A3j) 
will turn useful in the formulation of a general conjec- 
ture concerning the surface potential fa, see Appendix 
IBl We also note that for the systems under investigation 
here, the surface potential is quite high, and a very good 
approximation to (|A3|) is 
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APPENDIX A 

In order to explicitly relate the preferential coefficient 
r in ([2|) to the electric potential, we follow a procedure 
similar to that which leads to an analytical solution in 
the cell model, without added salt [24j. Implicit use will 
be made of the boundary conditions associated to (JTJ) . 
First, integrating Eq. ([1]), one gets 



(e - *+*-e*-*) r'dr', (Al) 



where the notation [-F(r')]^ = F(r) — F(a) has been in- 
troduced. Then, multiplying Eq. |T]) by r 2 tfi' and inte- 



We start by analyzing a 1:1 electrolyte, for which it 
has been shown [l9|, [2(| that the short distance behaviour 
reads 
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2A*log(y)-2*(0)] + O(nr) 4 (Bl) 



where ^ denotes the argument of the Euler Gamma func- 
tion *(x) = arg[r(ix)] [H, S3]- In |(BIj| . (j, denotes the 
smallest positive root of 



tan [2^1og(Ka/8) - 2*(/x)] 



2fi 



(B2) 



Expressions (|B1[) and (|B2[) require that £ exceeds a salt 
dependent threshold [denoted £ c below and given by Eq. 
(|C5p ] that is always smaller than 1 [l8j]. They thus always 
hold for £ > 1 and in particular encompass the interesting 
limiting case £ = 1, which is sufficient for our purposes. 
For large £, we have proposed in [l8| an approximation 
which amounts to linearizing the argument of the tangent 
in (|Blj) in the vicinity of — tt, and similarly linearizing ^ 
to first order: $(x) ~ — ir/2 — jx + C(x 3 ) where 7 is 
the Euler constant, close to 0.577. It turns out however 
that finding accurate expressions for exp(— z + 0o) ; which 
is useful for the computation of the preferential interac- 
tion coefficient, requires to include the first non-linear 
correction in the expansion of the tangent. After some 
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algebra, we find : 



diverge. This yields the leading order behaviour : 
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where the constant C = C 1:1 reads C 1:1 = 7 — log 8 ~ 
-1.502 and ^ (2) (!) = rf3 lnr(a;)/rfa; 3 | x=1 . From JB3]) and 
(|Bip where the sinus is expanded to third order, we ob- 
tain 



M 2 e-^4[(£-l) 2 + M 2 ] 
where Jl is given by 



logM +c -(*+£- 1)- 
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In writing (|B5D . we have introduced the change of vari- 
able /i = 2/i [23|. The reason is that similar changes for 
other electrolyte asymmetries allows to put the final re- 
sult in a "universal" (electrolyte independent) form, see 
below. A similar reason holds for introducing 24., here 
equal to 1, in the denominator of (|B5[) . 

The functional proximity between our expressions and 
those reported in [9j in the very same context is striking. 
We note however that our Jl (denoted j3 in Q) involves 
a different constant C. More importantly, the functional 
form of (|B1[) differs from that given in ■ The compari- 
son of the performances of our results with those of Q is 
addressed below, and is also discussed in the main text. 

Performing a similar analysis as above in the 1:2 case 
where z+ = 2 and z_ = 1, we obtain from the expressions 
derived in fijl: 



(ko) 



-z+<Po 



3[«-l) 



(B6) 



and similarly, in the 2:1 case (z+ = 1, z_ = 2): 

M 2 e-^° ~6[(£-l) 2 + £ 2 ]. (B7) 

In both cases, provided again that £ is not too low (see 
below) /i is given by (|B5[) [261 ] . with however a different 
numerical value for C [C 1:2 = 7 - (3 log 3) /2 - (log 2) /3 ~ 
-1.301 and C 2:1 = 7 - (31og3)/2 - log2 ~ -1.763]. 

The similarity of expressions (|B4[) , (|B6[) and (|B7|) leads 
to conjecture that this form holds for any z_:z+ elec- 
trolyte : 
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We then have to determine the prefactor A as a function 
of z + and z_. To this end, we make use of the exact 
relation (|A3[) [or equivalently (|A4|) ]. where in the limit of 
large £, the Ihs is finite while the two terms on the rhs 
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z + 

It then follows that A = 2(z + + z_)/z + so that our gen- 
eral expression (|B8|) takes the form: 
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This expression holds regardless of the approximation 
used for pt. If Eq. (|B5[) is used, then z + £ should not 
be too close to unity (see appendix [C] for more general 
results including the case z + £ = 1). 

In order to test the accuracy of (|B10|) in conjunction 
with (|B5p , we have solved numerically Eq. (TT|) for several 
values of na < 1 and electrolyte asymmetry and checked 
that for several different values of z + £ > 1, the quantity 



Q = -it 
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2{z+ + z_) 
log(Ka) + (z+£_ - I)" 1 
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is a constant C, which only depends on z + /z_ but not on 
salt and £ [it should be borne in mind that Eq. (|B5[) is a 
small na and large £ expansion, which becomes increas- 
ingly incorrect as na is increased and/or £ lowered]. This 
is quite a stringent test (since the two terms on the rhs of 
(jBlip are large and close] which requires high numerical 
accuracy. This is achieved following the procedure out- 
lined in (2?| • In doing so, we confirm the validity of (|B10j) 
and collect the values of C given in Table |TJ In the 1:1 
case, we predict that C — 7 — log 8 ~ —1.507, in excellent 
agreement with the numerical data of Figure [6l On the 
other hand, the prediction of Ref. 0] that Q reaches a 
constant close to -1.90 (shown by the horizontal dashed 
line in Fig[6]) is incorrect. Figure [6] shows that the quality 
of expression (|B10[) deteriorates when na increases, as ex- 
pected. It is noteworthy however that for na = 10 _1 , its 
accuracy is excellent whenever £ > 2. The inset of Fig. 
[6] shows the validity of (|B10[) for a 3:1 electrolyte. When 
z+t; is close to 1, Eq. (|B5[) becomes an irrelevant approx- 
imation to the solution of (|B2|) . and can therefore not be 
inserted into the general formula (|B10[) . This explains 
the large deviations between Q and the asymptotic value 
C observed in Fig. [S] for the lower values of £ reported. 
We come back to this point in Appendix [Cl 

The present results hold for z + £ > 1 + 0(1/ 1 logrea|). 
In this regime, our analysis shows that Eq. |B10p [with Jl 
given by (|B5|) ] is correct up to order 1/log (ko) for any 
(z_,z + ). On the other hand the results of [9( , valid in 
the 1:1 case, appear to be correct to order l/log 2 («;a). 
In addition, our expression for the surface potential may 
be generalized to a broader range of £ values and an ex- 
pression for the short distance dependence of the electric 
potential may also be provided. This is the purpose of 
appendix [Cj 
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FIG. 6: Plot of the quantity Q defined in (|B11[) versus line 
charge £ for a 1:1 electrolyte at ko = 10 -3 (continuous curve) 
and no. — 10~ 1 (dashed curve). The value reached at large £ is 
compared to the prediction of 9] Q — * e 7 + log 2 — 7 ~ —1.90 
(horizontal dashed-dotted line) whereas Eqs. (|B10|) and (|B5[) 
imply Q — > 7 — log 8 ~ —1.50, shown by the horizontal dotted 
line. The inset shows the same quantity for a 3:1 electrolyte 
at (ta = 10~ 5 [such a very low value is required to determine 
precisely the value of the asymptotic constant C, that can 
subsequently be used at experimentally relevant (higher) salt 
concentrations]. Here, we obtain Q — * —1.94 (dotted line) 
which is the value reported for C in Table [I] 



APPENDIX C 

In Appendix [B] the "universal" results valid for all 
(2+, z_) have been unveiled partly by a change of variable 
fi — > Jl from existing expressions [l8j. In light of these 
results, and of their accuracy (assessed in particular by 
the precision reached for the preferential interaction co- 
efficient), it is tempting to go further without invoking 
approximations of ()B2|) . or related expressions for other 
asymmetries than 1:1. Inspection of the results given in 
|18| for the 1:1, 1:2 and 2:1 cases lead, with again the 
help of (|A4[) . to the conjecture that 



jl y + zJ) 
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with 



tan [/ilog(«:a) + /iC] = 



(C2) 



We emphasize that (IC1[) . much as (jBll) . is a short dis- 
tance expansion and typically holds for nr < 1 (hence 
the requirement that na < 1). In appendix [Pi we give 
further analytical support for conjecture (|C1[) . A typical 
plot showing the accuracy of (|C1[) is provided in the main 
text (Fig. O. For nr < 0.1, the agreement with the ex- 
act result is better than 0.1%, and becomes progressively 
worse at higher distances (20% disagreement at nr = 1). 



From (jCip . it follows that the integrated charge q(r) 
in a cylinder of radius r [that is q(r) = — r0'(r)/2] reads 



z + q(r) 



-1 + /x tan 



/i log 
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where the so-called Manning radius [18j, l2a |29( is given 
by 



kR 



M 



exp — C 



2/1 



(C4) 



The Manning radius is a convenient measure of the coun- 
terion condensate thickness. It is the point r where not 
only z+q(r) = 1 but also where q(r) versus logr exhibits 
an inflection point [3(j. For high enough £, the logarith- 
mic dependence of 1/ju with salt [see (|B5I) ] is such that 
R M oc k- 1 / 2 . 

The two relations (|C1|) and (|C2[) encompass those given 
in Appendix [Bl and allow to investigate the regime z+£ c < 
and in particular the case z + £ = 1, the so-called 
Manning threshold |[. However, flUI]) and JCU) are not 
valid for £ < £ c , with 



1 



log Ka + C 



(C5) 



Note that ^ c < 1, since the constant C is negative and 
that salt should fulfill kcl < 1. For rea = 10~ 2 and 
= 1, we obtain £ c ~ 0.836. This is precisely the 
point where — T = 1 in the inset of Fig. [3] This inset 
also shows that the value of T resulting from the use of 
the solution of (|(J2[) is remarkably accurate. 

At this point, it seems useful to investigate the Man- 
ning threshold case z + £ = 1 (which corresponds to 
the onset of counterion condensation when na — > 
[1 [H, HI). It is readily seen that the solution of ([C2]) 
reads 



z+ {=i 



-tt/2 



log(Ka) + C ' 



(C6) 



which should be inserted in (|C1[) to obtain the potential 
profile, or in §5§ to get the interaction coefficient. 



APPENDIX D 

In this appendix we give further support for the con- 
jecture (jClj) which gives the short-distance expansion 
of the electric potential. Let us suppose initially that 
the charge is below the Manning threshold £ < £ c . 
It is straightforward to verify that Poisson-Boltzmann 
equation |T]) admits solutions which behave as <f>(r) — 
— 2A\n(nr) + InB + o(l) for nr <c 1. Injecting this ex- 
pansion into equation |T]) allows us to compute higher 
order terms. To study the regime beyond the Man- 
ning threshold, we compute all higher order terms of the 
form r 2n ( 1 + z + A ) (for a negatively charged macroion) and 
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r 2n{x-z-A) j-£ or a positively charged macroion), with n a 
positive integer. These terms turn out to present them- 
selves as the series expansion of the logarithm, thus re- 
summing them we obtain 



-2A ln(/tr 
2 

Z4 



hi 



1 



-InB 

z+B- z + ( Kr )2(i+*+A) 



(Dl) 



-A ln 

z_ 



1 



8(z+ + z_)(l + z+A) 2 

z_B z - (kt-) 2 ^ 2 -- 4 ) 
8(z + + z_)(l-z_A) 2 



The dots represent terms of order r 2 ™( 1 + 2 +- 4 )+ 2m ( 1 - z -- 4 ) 
with n and m two nonzero positive integers. When the 
Manning threshold is approached, z + A + 1 = for neg- 
atively charged macroion, the terms r 2n ( 1 + z + A ) (second 
line of Eq. (|D1[) ) become of order one, but the rest of the 
terms (third line of Eq. (|D1[) and dots) remain higher or- 
der: a change in the small distance behavior of <j) occurs. 
A similar situation is reached for 1 — z_A — which is 
the Manning threshold for a positively charged macroion. 

A and B in the previous equations are constants of in- 
tegration, which should be determined with the bound- 
ary conditions r<fr'(r) = 2£ at the polyion radius (r = a) 
and <fr — > for r — > oo. Thus to proceed further, we 
have to connect the long and the short distance behavior 
of <p. This connection problem has been only solved in 
the cases 1:1, 1:2 and 2:1 in Refs. [l^,[3l|. In particular, 
once A has been chosen (notice that for a = 0, A = — £), 
B should be one and only one function of A in order to 
satisfy cf> — > for r — > oo. The results from 0, [3l| show 
that 



B = 2 bA 7 ((l + A)/2) z (1 : 1) 
B 
B 



(D2) 



Z iA 2 ZA 1 (2(1 + A)/3) 7 ((1 + A)/3) (1 : 2) 

3 3A 2 2A 7 ((1 + 2A)/3) 7 ((2 + A)/3) (2 . x) 



where j(x) — T(x)/T(l — x). B turns out to have some 
interesting properties in the cases 1:1, 1:2 and 2:1, where 



its exact expression (|D2I) is known. Namely, at the Man- 
ning threshold 1 + z + A = 0, 



lim 



z+B~ 



A^-i/z+ 8(z + + z_)(l + z + A) 2 
Furthermore if we put 1 + z + A = ijl, and define 

z+B~ z + 



1 (D3) 



,2i*(£) _ 



8(z+ + z-)(l + z+A) 2 



(D4) 



then for /i G M, , 3 / (/i) 6 1 is a reaZ function of ju, with 
*(0) = 0. 

Let us now study the regime beyond the Manning 
threshold for a negatively charged macroion. From 
Eq. (|D1[) we can write 



,z+<t>(r)/2 



{ K r)- Z + A B Z +/ 2 f 1 - 



V ~ 8(z+ + z_)(l + z+A) 2 ; 



(D5) 

neglecting terms of higher order when z + A is close to —1. 

Let us conjecture that the properties of B as a function 
of A presented above hold in the general case z_ : z + . 
Then using the parameter Jl defined above we find after 
some simple algebra 



M V 2 ( Z + + z - 



sin [/xlog(«r) + ^(fJ.)] 
(D6) 



Recalling that \Jl\ <C 1 we can approximate ^(Jl) — jlC, 
where C = ^'(0). Replacing this approximation into (|D6|) 
and imposing the boundary condition a(f>'(a) = 2^ leads 
to (|C1|) and (|C2|) . Numerical values obtained for the 
constants C are reported in Table Q] for different charge 
asymmetries z_ : z + . The previous analysis shows that 
analytical predictions for C could be made if the con- 
nection problem is solved and the equivalent of expres- 
sions (|D2I) are found for the general case z_ : z + . 
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